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LONG-TERM  GOAL 

The  full  use  of  Lagrangian  observations  in  ocean  prediction  could  significantly  advance  the  Navy's 
ability  to  predict  ocean  conditions.  The  primary  goal  is  the  development  of  a  fully  operational, 
integrated  data  assimilation  scheme  will  afford  a  full  naval  predictive  capacity  in  fixed  ocean  regions 
which  can  be  comprehensively  surveyed  by  Lagrangian  measuring  devices.  A  key  part  of  this  aim  will 
be  to  develop  the  ability  to  assess  efficient  Lagrangian  deployment  systems.  This  will  be  based  on  the 
use  of  dynamical  systems  ideas  that  can  generate  optimal  strategies  for  deploying  Lagrangian 
observational  devices  and  their  associated  sensors.  The  optimal  deployment  coupled  with  an  effective 
Lagrangian  data  assimilation  strategy  will  form  the  basis  of  an  integrated  prediction  scheme  for  the 
ocean  that  can  feed  on  both  purely  Lagrangian  and  mixed  source  data. 

OBJECTIVES 

This  project  aims  to  develop  an  operational  technology  for  assimilating  Lagrangian  data.  This  new 
Lagrangian  data  assimilation  platfonn  is  expected  to  be  particularly  effective  in  ocean  regions  where 
coherent  structures  such  as  ocean  eddies  dominate  the  circulation.  The  focus  is  on:  1)  The  extension  of 
our  Lagrangian  data  assimilation  (LaDA)  approach  into  a  flexible  platform,  through  which  a  variety  of 
moving  instrument  measurements,  which  may  not  be  viewed  as  purely  Lagrangian  in  a  conventional 
sense,  can  be  integrated;  2)  The  design  of  observing  systems  that  take  full  advantage  of  all  moving 
instruments;  3)  The  formulation  of  automated  algorithms  for  optimal  deployment  strategies  of  the 
moving  instruments  so  as  to  maximize  the  information  content  of  the  observations;  and  4)  The 
incorporation  of  dynamical  systems  theory  to  enhance  our  predictive  skill,  in  particular  trough 
deciphering  coherent  structures  and  tracer  fields  associated  with  them. 

We  use  the  tenn  “Lagrangian  data”  in  the  broadest  sense,  including  the  positions  of  various  moving 
instruments  along  the  trajectories  as  well  as  time  sequence  of  the  tracer  fields.  Autonomous 
underwater  vehicles  (AUVs)  that  glide  or  maneuver  in  the  ocean  are  new  types  of  moving  instruments 
that  will  be  incorporated  into  our  extended  LaDA  platform  as  a  part  of  the  control  variables.  AVs  have 
come  to  be  a  stable  and  reliable  means  to  measure  water  properties.  The  platform  will  be  ideal  for 
designing  an  autonomous  ocean  sampling  network,  adaptive  observations,  and  optimal  deployment 
plans  of  such  moving  instruments.  The  extension  will  also  have  the  capability  to  handle  high- 
dimensionality  and  nonlinearity  of  the  operational  ocean  models.  An  example  of  the  structure  of  the 
tracer  fields  is  an  ocean  color  obtained  by  satellite. 
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A  comprehensive  LaDA  platform  will  be  built  upon  balanced  integration  of  data  assimilation 
techniques  and  dynamical  systems  concepts.  Recent  developments  in  the  theory  and  use  of  dynamical 
systems  concepts  for  Lagrangian  analysis  will  add  a  new  dimension  to  the  LaDA  platfonn. 
Incorporating  the  geometrical  approach  will  offer  a  much  needed  template  for  the  observing  system 
design  of  the  LaDA  platform.  In  return,  accurate  estimation  and  forecast  of  the  Eulerian  flow  field  by 
assimilating  such  data  will  further  strengthen  the  reliability  of  the  observing  system  with  moving 
instruments. 

APPROACH 

Our  approach  is  to  use  LaDA  as  a  basis  for  the  development  of  an  operational  technology  that 
accommodates  the  assimilation  of  data  from  a  variety  of  measurements  by  any  type  of  moving 
instruments  including  drifters,  floats,  and  AUVs.  Such  a  platform  is  being  developed  through  a 
hierarchy  of  models.  Fundamental  issues  are  addressed  using  idealized  model  flows.  Basic  tests  are 
being  performed  on  intennediate  model  flows,  the  double-gyre,  shallow-water  model  with  various 
wind  forcings,  and  progressing  to  the  realistic  general  circulation  model  for  operational  application  in 
the  Gulf  of  Mexico.  We  will  investigate  to  what  extent  Lagrangian  data  can  be  used  to  aid  in  the 
estimation  of  the  three-dimensional  flow  evolution.  We  will  improve  the  LaDA  method  to  better  deal 
with  the  chaotic  nature  of  the  Lagrangian  dynamics.  We  address  the  localization  and  inflation  issues 
that  will  help  simplify  estimates  of  error  covariance  and  develop  a  Lagrangian  Ensemble  Kalman  Filter 
(EnKF). 

Ideas  from  dynamical  systems  theory  will  be  employed  at  every  stage  to  provide  a  solid,  mathematical 
foundation  while  casting  conceptual  simplicity  on  the  complex  platfonn.  This  project  is  a  joint  project 
by  the  two  co-PIs,  K.  Ide  (UCLA)  and  C.K.R.T.  Jones  (UNC-CH).  G.  Vernieres  is  a  postdoctoral 
fellow  who  works  on  the  LaDA  platform  for  the  Gulf  of  Mexico.  H.  Salman  is  a  postdoctoral  fellow 
who  works  on  the  applications  of  the  dynamical  systems  theory  to  LaDA.  Liyan  Liu  is  a  graduate 
student  who  researches  the  vertical  propagation  of  the  information  by  the  LaDA.  Amit  Apte  is  a 
postdoctoral  fellow  working  jointly  with  Andrew  Stuart  (University  of  Warwick,  UK)  on  applying 
Langevin  sampling  techniques  to  Lagrangian  data  assimilation. 

WORK  COMPLETED 

Progress  is  being  made  in  advancing  Lagrangian  data  assimilation  towards  operational  use.  This  is 
being  achieved  through  a  coordinated  effort  to  test  and  adjust  the  LaDA  scheme  together  with  designs 
and  tests  of  optimal  deployment  strategies. 

Issues  essential  to  the  smooth  functioning  of  LaDA  have  been  investigated  during  the  past  year.  Due  to 
the  very  nature  of  Lagrangian  motion,  the  dynamics  is  usually  chaotic  at  least  in  many  sub-regions, 
filter  divergence  can  ensue  in  a  manner  that  is  a  challenge  to  anticipate.  The  underlying  cause  is  related 
to  the  distinguished  hyperbolic  trajectories  that  govern  the  flow  template  and  make  the  local  tracers 
separate  at  an  exponential  rate.  A  sophisticated  quality-control  method  is  being  developed  by 
combining  the  concepts  of  the  dynamical  systems  theory  with  an  expectation-maximization  (EM) 
algorithm. 

Much  Lagrangian  data  is  collected  with  surface  drifters.  A  natural  question  arises  as  to  how  effectively 
this  data  can  be  used  in  estimating  the  full  three-dimensional  ocean  flow.  We  have  implemented  LaDA 


in  a  two-layer  point  vortex  model.  The  potential  for  the  vertical  propagation  of  the  observed 
information  is  examined  through  testing  its  efficacy  in  parameter  estimation. 

Towards  the  development  of  an  operational  strategy  for  assimilating  Lagrangian  data,  the  strategy  for 
optimal  drifter  deployment  is  investigated  further  based  on  the  flow  template  regulated  by  the  invariant 
manifolds  using  the  EnKF-based  LaDA.  Currently  the  scheme  is  applied  to  a  single-layer  shallow- 
water  ocean  model.  We  have  discovered  that  a  mixed  strategy  of  placement,  near  hyperbolic  and 
elliptic  points  in  the  flow,  is  the  most  effective  when  the  observational  data  is  assimilated  through 
LaDA. 

We  have  implemented  a  multi-layer  reduced  gravity  model  of  the  Gulf  of  Mexico  (GOM.)  This 
modeling  setup  has  been  shown  to  be  the  simplest  representation  of  the  GOM  that  simulates  the 
shedding  of  eddies,  known  as  the  key  dynamic  effect  in  the  GOM.  It  is  a  first  step  toward  a  detailed 
representation  of  the  physics  of  the  GOM.  The  horizontal  discretization  uses  a  curvilinear  grid.  It  has  a 
horizontal  resolution  of  about  5  km  in  the  region  of  the  loop  current.  We  impose  a  sinusoidal  shaped 
current  at  the  Yucatan  channel  and  the  strait  of  Florida  with  respective  widths  of  160  km  and  150  km. 
In  order  to  conserve  the  total  mass  in  the  Gulf,  the  inflow  transport  matches  the  outflow  transport.  This 
configuration  allows  us  to  reduce  the  number  of  degrees  of  freedom  to  a  minimum  without  sacrificing 
the  resolution  near  the  region  of  interest.  The  structured  curvilinear  grid  of  the  horizontal  discrete 
domain  was  created  using  elliptic  grid  generation  techniques.  The  boundary  of  the  grid,  illustrated  in 
Figure  1,  follows  the  coastline. 


Figure  1.  Surface  Lagrangian  and  Eulerian  fields  of  the  Gulf  of  Mexico  by  the 
multi-layer  reduced  gravity  model  [Colors  represent  the  sea  surface  height 
(Eulerian  field).  The  black  lines  are  625  drifter  trajectories  over  eight  days 
(Lagrangian  field).  The  green  and  red  dots  represent  the  departure  and  arrival 

points,  respectively.] 

For  the  assimilation  of  the  tracer  field,  we  have  added  a  functionality  to  compute  the  passive  tracer 
field  as  represented  in  Figure  2. 
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Figure  2.  Lagrangian  tracer  fields  of  the  Gulf  of  Mexico  by  the  multi-layer  reduced 
gravity  model  [Colors  represent  the  normalized  concentration  of  the  tracer.  The 
dashed  black  lines  represent  the  sea  surface  height.] 


RESULTS 

By  combining  the  EM-algorithm,  the  tracer  control  scheme  can  detect  the  hyperbolic  effect  more 
reliably  and  promises  to  lead  to  a  robust  LaDA  method. 

The  LaDA  method  has  been  shown  to  afford  a  successful  estimation  of  the  highly-idealized  multi-layer 
point  vortex  model  within  the  Rossby  radius  of  deformation.  The  barotropic  mode  plays  a  key  role  as 
expected.  Testing  on  a  two-layer  model  is  natural  as  for  the  multi-layer  model  configuration  of  Gulf  of 
Mexico,  the  current  is  confined  to  the  two  upper  most  layers  with  the  ratio  of  transports  between  the 
two  upper  layers  being  kept  constant.  The  upper  most  layer  carries  80%  of  the  total  volume  transport. 

For  the  development  of  theoretical  deployment  strategies,  the  problem  is  made  complicated  by  the  fact 
that  the  positions  of  Lagrangian  data  are  continuously  changing  with  time  by  the  underlying  flow  field. 


In  order  to  tackle  this  problem,  an  understanding  of  how  the  underlying  flow  field  governs  the  motion 
of  Lagrangian  data  is  required.  This  problem  is  most  naturally  addressed  using  recent  ideas  borrowed 
from  dynamical  systems  theory  to  extract  Lagrangian  coherent  structures.  These  structures  have  been 
studied  extensively  for  geophysical  flows  in  recent  years  and  are  understood  to  orchestrate  the 
evolution  and  notion  of  material  particles. 

Figure  3  (left)  illustrates  the  locations  of  several  key  structures  that  were  identified  for  our  double  gyre 
flow  experiments.  Also  shown  are  the  deployment  sites  of  sets  of  drifters  in  relation  to  these  various 
structures.  Figure  3  (right)  illustrates  the  convergence  of  the  data  assimilation  method  for  several 
different  sites.  For  reference  purposes,  a  case  with  no  data  assimilation  is  included.  The  figure 
compares  four  different  cases;  unifonnly  released  drifters,  drifters  released  at  saddle  trajectories, 
drifters  released  within  coherent  vortices,  and  finally  a  combination  of  saddles  and  centers.  In  all  the 
cases  considered,  nine  drifters  were  used  for  the  data  assimilation.  The  results  clearly  show  that  the 
uniform  launch  produces  the  optimum  convergence  although  at  the  expense  of  having  to  release  the 
drifters  in  regions  remote  from  one  another.  In  stark  contrast,  targeting  Lagrangian  saddle  points 
appears  to  produce  equally  accurate  results  at  long  times  but  with  only  three  different  release  sites  for 
the  nine  drifters.  This  occurs  because  Lagrangian  saddle  points  can  disperse  drifters  over  a  wide  region 
providing  effective  sampling  of  the  entire  flow  field.  Therefore,  by  exploiting  the  internal  dynamics  of 
the  governing  flow  field  through  the  targeting  of  Lagrangian  coherent  structures,  we  are  able  to 
optimize  the  performance  of  our  method  through  the  optimal  deployment  of  drifters. 


day  =  60  0 


X  io*  days 


Figure  3.  Results  for  directed  drifter  deployment  strategy.  Left  pane:  Lagrangian 
coherent  structures  superimposed  on  the  height  field  at  day  60.  Right  pane:  Error 

in  the  height  field.] 

The  LaDA-EnKF  works  effectively  for  ocean  data  assimilation  using  a  shallow-water  model  (Figure  2) 
by  applying  the  localization  radius.  An  advantage  of  this  method  over  other  existing  methods  is  that 
the  assimilation  time  interval  can  be  as  long  as  the  Lagrangian  correlation  time  scale. 

Judicious  computation  of  the  error  correlation  and  innovative  construction  method  for  error  correlation 
are  crucial  in  the  coastal  data  assimilation  system  using  the  variational  approach. 


IMPACT/APPLICATION 


Much  data  for  the  ocean  is  Lagrangian  in  nature  and  our  techniques  promise  its  effective  use  in 
predictive  models.  We  are  developing  a  refined  scheme  to  incorporate  this  data  into  models  of  specific 
ocean  regions.  The  approach  is  adaptable  to  the  incorporation  of  all  data  comng  from  moving 
instruments,  including  floats,  AUVs  etc. 

As  the  scheme  is  developed  to  the  point  of  efficacy  in  real  time,  its  use  in  naval  operations  will  have 
signficant  application.  Its  impact  will  be  further  enhanced  by  the  optimal  dployment  strategies  being 
uncovered  through  the  application  of  dynamical  systems  ideas. 

The  development  of  LaDA  in  a  realistic  ocean  model  of  GOM  is  being  carried  out  on  a  model  that  is 
designed  to  offer  a  smooth  transfer  to  operational  naval  models. 
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